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ABSTRACT 

We present the results of a search for variability in and near the core of the metal-rich, 
obscured globular cluster Terzan 5, using NICMOS on HST. This extreme cluster has approx- 
imately solar metallicity and a central density that places it in the upper few percent of all 
clusters. It is estimated to have the highest interaction rate of any galactic globular cluster. The 
large extinction towards Terzan 5 and the severe stellar crowding near the cluster center present 
a substantial observational challenge. Using time series analysis we discovered two variable stars 
in this cluster. The first is a RRab Lyrae variable with a period of ~0.61 days, a longer period 
than that of field stars with similar high metallicities. This period is, however, shorter than the 
average periods of RR Lyraes found in the metal-rich globular clusters NGC 6441, NGC 6388 
and 47 Tuc. The second variable is a blue star with a 7-hour period sinusoidal variation and 
a likely orbital period of 14 hours. This star is probably an eclipsing blue straggler, or (less 
likely) the infrared counterpart to the low mass X-ray binary known in Terzan 5. Due to the 
extreme crowding and overlapping Airy profile of the IR PSF, we fall short of our original goal 
of detecting CVs via Pa emission and detecting variable infrared emission from the location of 
the binary MSP in Terzan 5. 



1. Introduction 

The study of variable stars in globular clusters 
has a long, rich history. Studies of RR Lyrae 
variables, long period variables and Population 
II cepheids have addressed fundamental questions 
about stellar evolution, globular cluster evolution 
and the distance scale. More recently, observa- 
tions of variable stars such as millisecond pulsars 
(MSPs), cataclysmic variables (CVs) and eclipsing 
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binaries have been motivated by studies of stellar 
interactions and globular cluster dynamics. For 
example, during stellar encounters binaries are ex- 
pected to become more tightly bound and give the 
energy difference up to passing stars, a process 
that can support a cluster against core collapse. 

Observations with HST have been crucial in 
overcoming the significant crowding present in 
such dense stellar environments. For example, 
HST observations of the center of the massive 
globular cluster 47 Tucanae have discovered two 
likely CVs (Paresce and de Marchi 1994 and Shara 
et al. 1996), nine eclipsing binaries (Edmonds et 
al. 1996) and 6 pulsating blue stragglers (BSs) 
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(Gilliland et al. 1998). Observations of the 
collapsed-core cluster NGC 6397 have discovered 
5 CVs (Cool et al. 1995, Grindlay et al. 1995; 
Edmonds et al. 1999; Grindlay 1999) and 3 likely 
low-mass WDs (Cool et al. 1998 and Edmonds et 
al. 1999) that may also be MSPs. 

HST has also been used to obtain the first 
high quality images and color magnitude diagrams 
(CMDs) of NGC 6441 and NGC 6388, two par- 
tially obscured clusters with surprisingly rich blue 
horizontal branches (HBs) when their high metal- 
licity is taken into account (Rich et al. 1997). 
Most metal-rich clusters have red, stubby HBs, in 
agreement with standard theories of stellar evo- 
lution. As metallicity decreases HBs usually be- 
come increasingly blue, but some clusters, such as 
NGC 6441 and NGC 6388, contradict this pattern 
and a second parameter is assumed to affect HB 
morphology. The cause of this "2nd parameter ef- 
fect" is a subject of vigorous debate, but some in- 
triguing clues have appeared. Both NGC 6441 and 
NGC 6388 appear to have HBs that are brighter 
than the HBs in metal poor clusters. Variability 
searches in these clusters have subsequently shown 
that the RR Lyrae stars have significantly longer 
periods than expected for metal-rich stars (Layden 
et al. 1999 and Pritzl et al. 2000), as does the RR 
Lyrae V9 in 47 Tuc (Carney et al. 1993). This ap- 
parently new subclass of long period, metal-rich 
RR Lyraes may show the lack of universality of the 
luminosity-metallicity relation for RR Lyraes and 
may require the establishment of a third Ooster- 
hoff group (Pritzl et al. 2000). 

The globular cluster Terzan 5, the subject of 
this paper, represents an excellent laboratory for 
the studies introduced above. Because its metal- 
licity is close to solar (Ortolani, Barbuy and Bica 
1996; Cohn ct al. 2000) it is ideal for studies 
of HB morphology and RR Lyrae behavior in a 
metal-rich globular cluster. It also has the highest 
predicted rate of stellar collisions of any galactic 
globular cluster (Verbunt & Hut 1987) making it 
a laboratory for stellar dynamics. It is already 
known to contain a LMXB (Johnston, Verbunt 
and Hasingcr 1995), a binary MSP (Lyne et al. 
1990; Nice & Thorsctt 1992), one other pulsar and 
a large number of unresolved pulsars (Fruchter & 
Goss 2000). Unfortunately Terzan 5 is highly red- 
dened, with E(B-V) = 2.49 (Ortolani ct al. 1996) 
and because of its large mass the crowding near 



cluster center, where the binaries are expected to 
be concentrated, is a formidable problem. 

Here, we describe an attempt to overcome the 
high reddening and crowding in Terzan 5 by ob- 
serving it in the IR with NICMOS on HST. Our 
aims were (1) to search the center of the cluster 
for CVs and the IR counterpart of the LMXB us- 
ing wide (F187W) and narrow-band (F187N) NIC- 
MOS filters centered on Paschen-a (Pa), and (2) 
to search for variability from the binary MSP. The 
general photometric results are described in Cohn 
ct al. (2000). In this paper we describe a detailed 
search for variability in two fields, one covering 
the center of the cluster and another containing 
the binary MSP. Although we fall short of detect- 
ing CVs because of crowding and no evidence is 
found for variability of any of the objects lying in 
the error circle of the MSP, we do detect an in- 
teresting RR Lyrae variable and a second variable 
that is likely to be an eclipsing BS, or (less likely) 
the IR counterpart of the LMXB. 

2. Observations and Analysis 
2.1. Data Reduction 

The Terzan 5 NICMOS observations described 
here were obtained on June 2nd and 3rd, 1998. 
Two pointings of 4 orbits each were obtained. 
Pointing (1) directed the NIC2 camera (19.2" x 
19.2" field of view) at the center of the clus- 
ter and pointing (2) directed the NIC1 camera 
(11.5" x 11.5" FOV) at the eclipsing MSP. Point- 
ing (1) was adopted for ^75% of each orbit and 
the telescope then slewed to pointing (2) for the 
remainder of each orbit. All three NICMOS cam- 
eras were used in parallel for these observations, 
but only the NIC2 data for pointing (1) and the 
NIC1 data for pointing (2) are of interest and are 
described here. 

Within pointing (1) a dither pattern was 
adopted, with different pointing positions for each 
of the 8 orbits. The pattern consisted of 8 out of 
9 points of a 6 pixel by 6 pixel square dither, with 
3 pixel steps. Superimposed on this box pattern 
were sub-pixel offsets in 1/8 pixel steps. Both of 
these factors were designed to improve the quality 
of the photometry, at the expense of the quality of 
the time series. The box pattern corrects the pho- 
tometry for image defects such as bad pixels and 
flat field errors and the sub-pixel offsets improve 
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the effective sampling of the images, for better res- 
olution. For the MSP pointing, no dithering was 
used to maximize the time series quality, but a 
different field position was used for the second set 
of 4 orbits to guard against accidental positioning 
on a bad pixel. 

For the central NIC2 observations, F110W, 
F187W and F187N were used for each of the 8 
orbits. The F187W and F187N filters are the in- 
frared analogues of R and Ha respectively and the 
F110W filter is similar to J. The latter was in- 
cluded to give color information for the observed 
stars. For the NIC1/MSP field only F140W im- 
ages were obtained. This filter was chosen because 
it gave the highest S/N and the greatest chance 
of detecting variability from the MSP, in partic- 
ular evidence of orbital variability from the sec- 
ondary, as in the system PSR1957+20 (Djorgov- 
ski & Evans 1988; Callanan, van Paradijs & Ren- 
gelink 1995). For the remainder of this section we 
will discuss only the central NIC2 observations. 

The NICMOS data were obtained in MULTI- 
ACCUM mode, which provides intermediate, non- 
destructive readouts during a single observation. 
For the F110W and F187W filters the 'STEP32' 
sequence was used, with 17 or 18 samples. The 
readout times were 0.3s, 0.6s, 1.0s, 2.0s, 4.0s, 8.0s, 
16.0s, 32.0s, 64.0s, etc (the remaining times in 
32.0s steps). For the F187N filter the 'STEP128' 
sequence was used with the same readout times 
as above except extending up to 128s readouts 
(with 22 samples). Standard processing was ap- 
plied via the STScI pipeline. Here, a problem 
was discovered with the processed F187N images, 
with a large number of stars found to contain pix- 
els with unusually small signal levels. With guid- 
ance from Mark Dickinson, Torstcn Boecker and 
Howard Bushouse at STScI, it was discovered that 
the saturation thresholds for the F187N process- 
ing were set incorrectly. Lowering the saturation 
thresholds by 5-10% and reprocessing using CAL- 
NICA repaired almost 100% of these faulty pixels. 

Customized software was used to convert the 
accumulated data into individual exposures suit- 
able for time series analysis. The data quality files 
were used to remove CRs and to avoid using satu- 
rated pixels. After blinking the individual images, 
pixels affected by 'grot' (probably small flecks of 
paint) were flagged as bad pixels and the bad pixel 
list was updated. Single bad pixels were corrected 



or 'cleaned' by adopting the average of the 4 ad- 
jacent pixels, but no attempt was made to correct 
groups of two or more adjacent bad pixels. 

2.2. Time Series Analysis 

Our star list for differential photometry was 
based on the work of Cohn et al. (2000) who av- 
eraged and combined the 8 separate images for 
each filter into single, over-sampled images using 
the STSDAS drizzle routine. The F110W filter, 
which gave the deepest, best resolved final image 
was used to calculate the master star list used here 
(containing 3091 stars). 

Two standard photometry techniques were used 
to produce time series: 

(1) Aperture photometry: A plate solution de- 
termined from a sample of bright stars was ap- 
plied at each time, to calculate the coordinate 
shifts to be applied to the master coordinate list. 
Aperture photometry was then carried out using 
DAOPHOT under IRAF, with an aperture ra- 
dius of 2.5 pixels, roughly matching the radius 
of the first Airy minima for the F187 data. To 
improve the signal-to-noisc (S/N) ratio of the in- 
tensity measurements and minimize the contribu- 
tion from near neighbors, the signal in each pixel 
in each aperture was weighted by a simple model 
of the PSF. Sky values were determined from the 
mode of the counts in an annulus around each star. 

(2) Point spread function (PSF) fitting: The 
extreme crowding of the central Terzan 5 field, 
especially the large number of bright giants, the 
complicated nature and large radial extent of the 
PSF and the marginal sampling of the individual 
images caused problems with PSF-fitting for the 
faintest stars. In particular, many faint objects 
in the master star list were not detected in the 
F110W images (the most undersampled data) us- 
ing ALLSTAR. Therefore, we experimented with 
a smaller sized list of 1702 stars containing ob- 
jects found in PSF-fitting of both the F110W and 
F187W images (Cohn ct al. 2000). 

We examined the counts for each star as a func- 
tion of time for the longest readout times (32s or 
128s), and found a high correlation between the 
time series of different stars, with a linear falloff 
in the counts as a function of time during an expo- 
sure, caused, we believe, by 'persistence' (Skinner 
et al. 1998). To remove this artifact we calculated 
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the average time series for bright non-saturated 
stars and normalized each time series by this av- 
erage function. Zero-points based on the photom- 
etry of Cohn ct al. (2000) were used. 

We then compared the two photometry tech- 
niques using the root-mean-square (rms) values 
calculated for each star. Only the exposures with 
the longest sample times (either 32s or 128s) were 
included in this calculation. Aperture photometry 
gave significantly better results than PSF fitting 
for the F110W data and comparable or slightly 
better results for the F187W and F187N data. 
Since the aperture photometry results were supe- 
rior we will quote them unless otherwise noted. 

Figure 1 shows plots of rms versus F110W, 
F187W and F187N magnitudes. The solid line 
shows the theoretical rms based on the limits 
from photon noise and readout noise (Kjeldsen & 
Frandsen 1992). We have plotted the average rms 
limit in 0.5 mag wide bins in each case. At faint 
magnitudes our measured rms values are only fac- 
tors of 2-3 above the Poisson limit, a reasonable 
result given the large amount of dithering adopted, 
and other systematic effects. 

Because the images are crowded we have ex- 
amined the degree of overlap between stars. We 
define the 'contamination' level (CL) for any star 
to be the amount of light from other stars falling 
within its aperture (for example, CL = implies 
a star is completely isolated and CL = 1 implies 
equal contribution from the star and neighbors). 
The CLs for the F110W filter are shown in Figure 

2. The median CL is 1.70, showing the high level 
of crowding in the data. 

To search for variables we concentrated on time 
series analysis. We placed the stars in order of 
decreasing brightness, divided the stars into bins 
containing 50 stars and searched for stars with un- 
usually large rms values in each bin, using an it- 
erative procedure. Candidate variables were iden- 
tified as having unusually large (> 3cr) rms values 
in all three of the filters. 

3. Results 

3.1. Detected variables 

Thirteen candidate variables were discovered 
using the above method. All but three of these 
stars have CL > 2.0 in the F110W image. The 



second most isolated candidate variable (at 194.05, 
197.78 in orbit #1 coordinates) has CL=0.37 and 
is clearly a good candidate variable. This star will 
be referred to as VI and is located at a(2000) = 
17 h 48 m 05 s . 14, 5(2000) = _24°46'38".2. Figure 1 
shows VI along with stars lying within 12 pixels 
of VI ('a'-'e'). Star 'a' is a faint star lying only 2.1 
pixels from VI. The overlapping PSF of the vari- 
able imposes an apparent but false signal for its 
fainter companion. None of the other stars have 
convincing variations. Further discussion of VI 
will be given in section 3.2. 

A deeper search was also carried out for 2er 
variations. This resulted in 90 variable candi- 
dates, with 13 stars having CL< 2.0, but this 
deeper study revealed only one new, believable 
variable, henceforth V2 (see section 3.3), lo- 
cated at a(2000) = 17 h 48 m 05 s .06, 5(2000) = 
-24°46'52".5. A finding chart (Figure 3) shows 
the F110W NIC2 image from the first orbit and 
the positions of VI and V2. The summed im- 
age from the first orbit is shown (rather than the 
dithered, oversampled image) to show the level of 
crowding in the single orbit images. 

To test the completeness and depth of our vari- 
able search, we carried out time series simulations 
of our F110W data, the least crowded images. We 
calculated aperture photometry time series for a 
grid of positions and added sinusoids with an am- 
plitude of 0.3 mag to these time series. We con- 
sidered a fake variable to be "detected" if the rms 
was more than twice the rms without adding in 
the sinusoid. This technique automatically takes 
crowding into account. We then tested detection 
efficiency as a function of radial distance and ab- 
solute magnitude of the fake variables (for periods 
between 0.5 hours and 15 hours, the detection ef- 
ficiency is constant to within a few percent). The 
My calibration was based on the assumption that 
the simulated stars belong to the main sequence or 
the giant branch, using stellar isochrones (see Sec- 
tion 3.2.1). The percentage of detected 0.3 mag 
simulations as a function of radial distance from 
the center of the cluster is shown in Figure 4. Note 
the dramatic improvement in detection efficiency 
as one moves away from the center of the clus- 
ter, but the lack of sensitivity to variables with 
My > 4 (just above the main sequence turn-off; 
MSTO). 
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Variable searches were also attempted using 
analysis of the Lomb-Scargle periodogram (Scar- 
gle 1982) of each time series. However, the 8 or- 
bits were spread out over ^1.3 days, and the noise 
from orbit to orbit is larger than the noise within 
a single orbit, because of the dithering. Therefore, 
many false signals with periods of a few hours were 
found in the power spectra and this method proved 
to be less efficient at finding candidate variables 
than working in the time domain. 

3.2. VI: RR Lyrae Variable 

We applied several extra tests to VI to confirm 
its variability. Because a bad pixel at (197,203) 
falls within the VI aperture for the 4th and 5th 
dither positions (1.2 and 1.5 pixels away respec- 
tively), we verified that there is no clear drop in 
the flux for these two dither positions. Besides 
these two dither positions the average offset of 
the bad pixel from VI is 3.9 pixels, well outside 
our aperture. As a second test, we placed a grid 
around VI and calculated time series for apertures 
centered on this grid. Variability was only found 
for apertures within a few pixels of VI. Finally, we 
verified that VI 's light curve from PSF-fitting was 
similar to that for aperture photometry, and that 
the quality of the fits was reasonable compared to 
other stars. 

Henceforth we consider VI to be a true variable. 
We now discuss the photometry and time series for 
this star. 

3.2.1. Photometry 

Figure 5 shows the F110W vs F110W-F187W 
CMD for the central NIC2 field, from Cohn et al. 
(2000). The labels represent the variable VI and 
its near neighbors ('a'-'e'). The variable is clearly 
a blue star, while the stars labeled with '*' are 
other possible blue stars, many of which are likely 
to be red giants scattered into the blue region of 
the CMD by errors. The solid line shows the 
stellar isochrones of Bergbusch and Vandcnbcrg 
(1992). We extrapolated the highest metallicity 
stellar isochrones ([Fe/H] = -0.47) to solar metal- 
licity, fed the T e // and log g values into Kurucz 
atmospheres and used SYNPHOT to estimate the 
NICMOS magnitudes of these Kurucz models, us- 
ing the reddening measured by Cohn et al. (2000). 
Besides the main sequence and giant branch we 



also show simple models of the horizontal branch 
(HB), assuming a large spread in effective temper- 
ature. Since the variable lies close to an extension 
of the HB, a RR Lyrae source for the variations 
is suggested. We also show a model BS sequence 
(guided by T e f f and log g values from Shara, Saffer 
and Livio (1997) for a bright BS in 47 Tuc). Our 
BS sequence extends to ^2.2 magnitudes above 
the MSTO in an unreddened V vs B - V CMD, 
about equal to the extent of the BS sequence in 
the center of 47 Tuc obtained with HST (Gilliland 
et al. 1998). 

Despite the very small FOV of the NIC2 cam- 
era there is some field star contamination because 
the line of sight towards Terzan 5 is close to the 
galactic center (£=3.81, B=1.67). We estimated 
the field star component using the ground-based 
observations of Ortolani et al. (1996). Using their 
SUSI (2.2' x 2.2' FOV) observations extending to 
about 1.5 magnitudes below the HB (see their Fig- 
ure 5) we defined field stars to have V — I < 3.2. 
Scaling to the NIC2 FOV we predict 4-5 bright 
field stars should be found in our field. This is 
in excellent agreement with Figure 5 showing 5-6 
likely field stars with F110W < 19.0 and F110W- 
F187W <-0.8. 

3.2.2. Time series 

Figure 6 shows the aperture photometry light 
curves of VI along with light curves of stars 'a'- 
'e'. The F187W and F187N light curves have 
been shifted vertically so that the apparent peaks 
of the light curves (in orbit 3) match up. Be- 
sides the overlapping star 'a', VI clearly stands 
out from its near neighbors. A distinctive element 
of the light curve of VI is its asymmetry, typi- 
cal behavior for RR Lyraes of type ab. Motivated 
by this, we carried out fits of our NICMOS light 
curves to two sets of IR light curves for RR Lyraes, 
if -band templates from Jones, Carney and Ful- 
bright (1996) and the JHK light curves of Cac- 
ciari, Clementini and Fernley (1992) for three field 
RRab Lyraes. We varied the period, magnitude 
zero-point and amplitude of the IR light curves 
(tabulated as functions of phase) to minimize the 
residual between our light curves and the tabu- 
lated light curves. Since the F187W and F187N fil- 
ters lie roughly halfway between the ground-based 
H and K passbands we used both H and K data 
to fit to our combined long- wavelength (F187W 
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and F187N) light curves. We used J data to fit to 
our F110W light curves. 

Jones ct al. (1996) give templates for (a) RRc 
variables, (b) RRab variables, A(B) < l.Omag, 
(c) RRabs, 1.2 < A(B) < 1.3, (d) RRabs, A(B) « 
1.5, and (e) RRabs, 1.6 < A(B) < 1.7. By exam- 
ining the residuals and comparing them with the 
expected scatter we can reject (a) at the ~ 3. So- 
level (the derived amplitude of 0.21 mag is also too 
high by about a factor of 2). The best fits were 
obtained for the (c), (d) and (e) templates, with 
little difference between these three. With (d) for 
example, we derived a period of 0.60 days and a 
peak-to-peak amplitude of 0.25 mag when fitting 
to the 1.87/xm filters. 

The Cacciari et al. (1992) light curves for W 
Tuc gave the best fit to our data. For J we de- 
rived a period of 0.63 days and an amplitude of 
0.44 mag when fitting to F110W. Using H and 
K time series when fitting to the 1.87/zm filters, 
we derived periods (amplitudes) of 0.61 days (0.24 
mag) and 0.61 days (0.25 mag) respectively. The 
ratio between the amplitudes derived for J and for 
K was only 15% larger than the ratio for the field 
systems, a reasonable agreement given the differ- 
ences between the HST and ground-based filters. 
The similarity between the derived periods shows 
that the error for the period determination is likely 
to be only 5%-10%. The rms of the residuals from 
these fits are 0.018 (Jones K), 0.021 (Cacciari J), 
0.017 (Cacciari H) and 0.016 (Cacciari K), consis- 
tent with the noise properties of stars in a similar 
magnitude range (see Figure 1). 

Figure 7 shows the results of these fits. Note the 
similarity between the H and if-band template 
light curves, the close correspondence between the 
F187W and F187N data and the rapid rise time 
for the F187N data near 7000s (matching well with 
the template light curves). 

3.3. A Search for Blue variables: V2 

Along with the general search for variables in 
the NIC2 images, we also paid special attention to 
stars with blue colors, in an attempt to find the in- 
frared counterpart to the LMXB known in Terzan 
5, and to search for other RR Lyrae variables and 
possible bright BS variables. 

To select a robust sample of stars for time se- 
ries analysis, only the stars with CL < 1.5 in the 



F187W image were examined in greater detail. To 
compare our photometry with that of Cohn ct 
al. (2000), we show in Figure 8(a) the F110W 
vs F110W-F187W CMD and in Figure 8(b) the 
F187W vs F187N-F187W CMD based on the flux 
values derived from our time series. The labeled 
stars from Figure 5 are shown, plus extra stars 
with apparently blue colors from our photometry 
(labeled with '+'). Note that: (i) the RR Lyrae 
variable is found in a similar part of the F110W 
vs F110W-F187W CMD as Figure 5, (ii) many of 
the blue stars labeled '*' from Figure 5 are missing 
from Figure 8 because they are too contaminated, 
and (iii) of the blue stars that do appear, most 
of them, reassuringly, have blue colors. Apparent 
Pa-bright objects in Figure 8(b) are caused by er- 
rors from near neighbors. 

Figure 9 shows the corresponding plots of rms 
versus magnitude. Using an uncontaminated sam- 
ple of stars causes a significant reduction in the 
number of false variables appearing in these plots. 
In particular the RR Lyrae variable and its near 
neighbor stand out much more clearly from the 
other stars. The star labeled with 'V2' also has 
an unusually large rms, especially in the F110W 
filter. This large rms is unlikely to be caused by 
near-neighbor effects, since this star is well isolated 
from other stars, having CL = 0.05, 0.24 and 0.22 
in F110W, F187W and F187N respectively. (The 
dithered images show a faint star about 2 pixels 
away from V2 that is too faint to be included in 
the master star list. However, this star makes < 
10% contribution to the counts in V2's aperture.) 
Also, there are no bad pixels or other obvious ar- 
tifacts near this star. Figure 10 shows plots of the 
time series for V2. 

The absolute magnitude of V2 is uncertain be- 
cause we lack optical information for this star and 
only have limited color information. Noting that 
V2 lies near the end of our BS sequence in Figure 
5, we assume that a 8000 K Kurucz atmosphere 
gives a reasonable approximation to the spectrum 
of V2, and thus estimate that My ~ 2.0. 

Power spectrum analysis gave periods and am- 
plitudes of (6.94 hr, 0.27 mag; F110W), (6.93 
hr, 0.33 mag; F187W) and (6.88 hr, 0.27 mag; 
F187N), assuming the variation is sinusoidal. By 
adapting the fitting technique used for the RR 
Lyrae variable (using a sinusoid for the template), 
we estimated the periods and amplitudes to be 
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(6.96 ± 0.20) hr and (0.27 ± 0.09) mag (F110W), 
(6.74 ± 0.25) hr and (0.40 ± 0.19) mag (F187W) 
and (6.85 ± 0.6) hr and (0.29 ± 0.27) mag (F187N) 
respectively, in good agreement with the power 
spectrum results quoted above. Using this fitting 
procedure, the rms of the residuals were found 
to be 0.045 mag, 0.062 mag and 0.122 mag for 
the three filters. These values agree well with the 
mean rms values (0.045 ± 0.010, 0.069 ± 0.026 and 
0.111 ± 0.032) determined for the bins containing 
V2 in Figure 9. While a slightly larger amplitude 
is suggested for F187W, the errors for this filter 
are quite large, and a single sinusoidal fit to all 
of the data is allowed by the errors. For example, 
adopting a 'compromise' solution for all 3 filters of 
period 6.98 hours and amplitude 0.282 mag gives 
rms residuals of only 0.049, 0.075 and 0.126 mag, 
comfortably within the ler limits quoted above. 

3.4. Millisecond pulsar 

3.4-1- Astrometry 

The position of the MSP in Tcrzan 5 is known 
to within 0.08" (Fruchter & Goss 2000) in the ra- 
dio frame of reference (which uses quasars), but 
positional information given in the NICMOS data 
files is accurate to only about 1" in absolute terms, 
and is limited by uncertainties in the locations of 
the instrument apertures relative to the Optical 
Telescope Assembly axis and the inherent uncer- 
tainty in the positions of the two guide stars used. 
For HST images, differential astrometry (measur- 
ing relative positions within an image) is much 
more accurate. We therefore carried out an astro- 
metric study of the Terzan 5 field to improve the 
absolute accuracy of the positions, and cut down 
the size of the error circle for the MSP. Our posi- 
tions should then only be limited by uncertainties 
in our astrometric solution and by systematic er- 
rors in the STScI Digitized Sky Survey (DSS). 

We used two steps to carry out this astrometric 
study: (i) We matched up the maximum possi- 
ble number of stars (9) between the DSS and the 
study of Ortolani et al. (1996), and used a linear 
regression in x and y to calculate an astromet- 
ric solution for the Ortolani stars, (ii) We then 
performed a plate solution match between the Or- 
tolani and NICMOS fields. A direct solution be- 
tween the DSS image and the NICMOS fields was 
not possible because the small FOV of NICMOS 



and poor resolution of the DSS images limited the 
number of matching stars. 

The position of the MSP is: 
a(2000) = 17 h 48 m 02 s .21, 5(2000) = -24°46'37".2 
(Fruchter & Goss 2000). Using the two steps de- 
scribed above, the nominal MSP position corre- 
sponds to pixel coordinates of (141.6, 180.5) for 
orbits 1-4, only 0.12" from the predicted position 
using the header information. The rms residuals 
after applying the fit described in step (i) above 
are 0.341" and 0.140" in the x and y directions 
respectively. The errors for step (ii) are much 
smaller. 

3-4-2. Photometry 

Nice and Thorsett (1992) estimate that if the 
MSP companion is a main sequence star it is likely 
that 0.09 M Q < M 2 < 0.11M©, right near the 
lower limiting mass for hydrogen burning. How- 
ever, by analogy with the eclipsing MSP 1957+20, 
the optical counterpart is expected to be much 
brighter (My ~ 10 — 11) than a normal main 
sequence star due to reprocessing of pulsar spin- 
down energy on the evaporating wind. A search 
for photometric variability from the MSP is nev- 
ertheless very challenging given the considerable 
distance of the Terzan 5 MSP and its high redden- 
ing. These problems were compounded by 3 extra 
problems: (1) we had initially planned to observe 
the cluster center and the MSP simultaneously us- 
ing the NIC1 and NIC2 cameras (respectively) in 
parallel, given the virtually identical value for the 
cluster center-MSP offset and the NIC1-NIC2 off- 
set. However, in phase 2 planning we found that 
we were unable to use the required telescope roll 
angle, and therefore decided to slew between the 
two targets, concentrating most of our time on 
cluster center and cutting our observing time on 
the MSP by a factor of ^4. This strategy effec- 
tively reduced the observed amplitude of any MSP 
variations that may be present. (2) The nominal 
position of the MSP lies close to a HB star and two 
bright giants (at relative distances of 0.55", 1.3" 
and 2.5"; see Figure 11). The extra background 
from the PSF halos of these bright stars compro- 
mised the MSP observations. (3) The NIC1 data 
is much more affected by transient artifacts, in- 
cluding 'persistence', than the central NIC2 data. 

Despite these limitations the MSP-field data 
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is less crowded than the central NIC2 data, en- 
abling variability studies of stars at and below 
the MSTO. Our photometric analysis was simi- 
lar to that carried out for the central NIC2 field 
except that we did not have the advantage of cross- 
checking candidate variables between different fil- 
ters. A plot of time series rms versus magnitude 
for the MSP field is shown in Figure 12. The HB 
stars identified from the V vs V — I CMD of Or- 
tolani et al. (1996) are labeled, as are stars (and 
in some cases PSF tendrils) lying within 1" of the 
nominal MSP position, and the two giants near 
the position of the MSP. Note that the stars near 
the MSP have slightly higher noise levels than av- 
erage because of the higher background near the 
bright giant stars. 

We carefully examined the results to search for 
evidence of variability, but no authentic variable 
candidates were found. We also carried out a 
second dedicated search for variability from the 
MSP by searching a grid of positions within a 3a 
(1.1") radius circle around the measured position 
of the MSP. Two criteria were used in this vari- 
ability search: excess rms in a small group of ad- 
jacent pixels and high linear correlation of the 8 
point time series with the predicted time series 
of the MSP (assuming a sinusoidal shape for the 
light curve and using the ephemcris of Nice and 
Thorsett 1992). No candidates for a variable MSP 
were found. 

To estimate the depth of our MSP variability 
search, we used our Kurucz atmosphcrc/synphot 
simulations described above. We simulated MSP 
variations using sinusoidal light curves (at the 
MSP period of 1.82 hr) with an amplitude of 1 mag 
and added these light curves to our measured light 
curves for the above grid. A MSP simulation was 
considered 'detected' if the linear correlation be- 
tween the predicted and simulated light curves was 
>0.5 and the rms of the simulated light curve was 
50% higher than the original rms. At My = 3.9 
(F140W=21.3) only 43% of the simulations were 
detected and at M v = 3.4 (F140W=20.0) 86% of 
the simulations were detected. To show the effect 
of the high sky levels near the expected position 
of the MSP we carried out another set of simula- 
tions for a grid centered on a point well away from 
the bright giants. Here, our variability search was 
able to go 1.5 mag deeper, but still falls well short 
of the desired levels to detect variability, assuming 



that the secondary has My ~ 10 — 11. 

4. Discussion 

4.1. RR Lyrae variable 

The identification of VI as a RR Lyrae variable 
and a cluster member is secure, based on the light 
curve, the position of the star in the CMD and the 
small chance of field star contamination. As with 
47 Tuc, only one RR Lyrae has been discovered in 
Terzan 5 and the bulk of the HB stars are red and 
lie near the giant branch (see Figure 5), as they 
do in most metal-rich globular clusters. However, 
Cohn ct al. (2000) report the possible detection 
of a small number of blue HB stars in Terzan 5. 
The high reddening and crowding in this cluster 
makes such detections much more difficult than in 
NGC 6388 and NGC 6441 where many blue HB 
stars have been found and there are multiple RR 
Lyrae variables. Rich et al. (1997) have explored 
two possibilities to explain the blue HB star popu- 
lations in NGC 6388 and NGC 6441, the effects of 
age and dynamics. At high metallicity an increase 
in age can populate a blue HB tail with stars at 
the low end of the mass distribution, while dynam- 
ical encounters in the cores of dense clusters can 
result in extra mass loss via tidal stripping. Rich 
et al. (1997) conclude that although neither effect 
alone appears sufficient to explain the cluster blue 
HBs, a combination of them may have an effect. 

Our discovery of a RR Lyrae variable in Terzan 
5 is therefore useful for the study of the influence 
of metallicity, age and interaction rate on the pro- 
duction of HB stars in the instability strip. For 
example, Terzan 5 has a higher metallicity than 
any of the above clusters which should result in 
fewer RR Lyraes, but it also has a higher inter- 
action rate which may result in a higher number 
of RR Lyraes (see Rich et al. 1997). It is possi- 
ble that these two factors combine. For example, 
maybe a large number of stars in Terzan 5 lose 
mass via tidal stripping but the cluster metallicity 
is high enough that only the stars with the most 
extreme mass loss show up as RR Lyraes. 

Sweigart and Catelan (1998) note that the HBs 
in NGC 6388 and NGC 6441 have a pronounced 
upward slope with decreasing B — V , a result 
which is not predicted by canonical HB mod- 
els and which cannot be reproduced by either a 
greater cluster age or enhanced mass loss along the 
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red giant branch. They explore three possibilities 
to explain these sloped HBs: (1) a high helium 
abundance, (2) a rotation scenario in which inter- 
nal rotation during the RGB phase increases the 
HB core mass, and (3) a helium-mixing scenario 
in which deep mixing on the RGB enhances the 
the envelope helium abundance. 

One prediction of the high-Y scenario is that 
the number ratio R of HB stars to RGB stars 
brighter than the mean RR Lyrae luminosity 
should be large. We calculated R for Terzan 5 
by defining the HB as those stars having 15.6 
< F110W < 16.1 and correcting for the large 
number of RGB stars found in this region of the 
CMD by averaging and subtracting the number of 
stars present in adjacent magnitude bins (15.1 < 
F110W < 15.6 and 16.1 < F110W < 16.6). Our 
value for R is 1.03 ± 0.11, where the error is based 
on Poisson statistics alone. Although systematics 
will increase the error, this value is clearly in- 
consistent with the high-Y scenarios presented in 
Sweigart and Catelan (1998), where R = 3 — 4. 
Layden et al. (1999) came to a similar conclu- 
sion for their study of NGC 6441. Sweigart and 
Catelan (1998) suggest a variety of other observa- 
tional tests for discriminating among the various 
scenarios, all requiring further observations. 

Figure 5 shows that VI may be somewhat 
brighter (~ 0.3 mag) than expected for a blue 
HB star, consistent with the sloped HB slopes 
mentioned above, although uncertainties in the 
large reddening for this cluster and the limitations 
of our simple HB model make detailed compar- 
isons difficult. As Sweigart and Catelan (1998) 
point out, unusually bright RR Lyraes should have 
unusually long periods, as observed for the RR 
Lyraes in NGC 6388, NGC 6441 and 47 Tuc. This 
self-consistency between the RR Lyrae luminosi- 
ties and periods is reassuring. 

The obvious question arises: is the period of 
VI unusually long? First, the metallicity has to 
be considered. VI is likely to be the most metal- 
rich RR Lyrae yet found in a globular cluster, since 
Terzan 5 is thought to have roughly solar metallic- 
ity (Ortolani et al. 1996) or even higher (Cohn et 
al. 2000 estimate [Fe/H]=0.25). NGC 6388, NGC 
6441 and 47 Tuc have [Fe/H] values of -0.60, -0.53 
and -0.76 respectively (Harris 1999) 2 . It may even 

2 The Catalog of Parameters for Milky Way Globu- 



be the highest metallicity RR Lyrae ever found. 

Examination of the metallicity-period relation- 
ship for field RR Lyraes (Layden 1995) shows that 
VI has a significantly longer period than field stars 
with similar metallicity. The period-metallicity 
distribution of Pritzl et al. (2000) shows that VI 
lies well outside the two Oosterhoff classes, with a 
longer than expected period given the metallicity 
of Terzan 5. Therefore VI joins the rapidly grow- 
ing subclass (suggested by Layden et al. 1999) 
of long period RR Lyraes in metal rich globular 
clusters. 

Also plotted on Figure 1 of Pritzl et al. (2000) 
are the periods and metallicities of the RR Lyraes 
in NGC 6388 and NGC 6441. These two clus- 
ters are clearly separate from the two Oosterhoff 
classes, with higher metallicity and longer periods. 
Interestingly, while the average period for these 
two clusters is longer than all of the Ool and OoII 
globulars, VI has a period lying almost exactly 
between the 2 Oosterhoff groups. Further, a line 
drawn between VI and the average of NGC 6388 
and NGC 6441 in the period-metallicity plot lies 
almost parallel to the line between the 2 Ooster- 
hoff classes, suggesting two period-metallicity re- 
lations with similar slope may exist. 

This interesting but speculative result clearly 
requires verification by observing a larger sample 
of RR Lyraes in clusters with solar metallicity or 
above. It could be that VI has an unusually short 
period compared to other (undetected) RR Lyraes 
in this cluster. For example, the periods in NGC 
6441 alone range from ^0.6 days to <~ 0.9 days, 
with a reasonably flat distribution between these 
two extremes. 

4.2. The faint variable V2: LMXB or 
eclipsing BS? 

To summarize our results, V2 is blue, relatively 
faint and variable, with a near sinusoidal light 
curve, and a period of ^7 hr. We consider two 
possibilities for this variable, that it is a LMXB or 
that it is an eclipsing BS. 

First, we consider the LMXB hypothesis. In 
this case the variability could be caused by at 
least three different phenomena: (1) Ellipsoidal 

lar Clusters, compiled by William E. Harris, McMas- 
ter University, revised June 22, 1999 is available at 
http://physun.mcmaster.ca/~harris/mwgc.dat 
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variations from distortion of the secondary. These 
variations typically have a distinctive double-lobed 
pattern, with a difference in depth between suc- 
cessive minima of ^0.05-0.10 mag and an ex- 
pected difference in amplitude between F110W 
and F187W of ^0.05 mag, using reasonable as- 
sumptions regarding mass ratios, temperatures 
and inclinations (van Hamme, private communica- 
tion). These amplitude differences are not clearly 
detected in the data, though the relatively large 
error for the F187W amplitude means we cannot 
rule out this possibility. 

(2) Another possible explanation of variability, 
under the LMXB hypothesis, is the effects of heat- 
ing of the secondary by the disk and neutron star. 
Here, the expected light curve is also double-lobed, 
but with the positions of the deep and shallow 
minima reversed (see van Paradijs & McClintock 
1995). A large contribution from heating would be 
obvious in the light curve because of the significant 
expected difference in eclipse depth, but a small 
contribution would plausibly cancel out the differ- 
ence between the depths of adjacent minima, when 
combined with an ellipsoidal variation. However, 
this requires some fine-tuning of the binary pa- 
rameters and a difference between the light curve 
amplitudes remains (Van Hamme, private commu- 
nication). Note that if either case (1) or (2) applies 
to V2 then the period would be twice the 7-hour si- 
nusoidal variability period, ie ~14 hours. It would 
thus be similar to the field LMXB, and recurrent 
transient, Aql X-l (Chevalier & Ilovaisky 1991), an 
attractive possibility given the transient, or highly 
variable nature of the Terzan 5 LMXB. 

(3) Flickering from the disk may also be ex- 
pected to show up in the light curve of a LMXB, 
as it docs in field systems. If V2 is a LMXB, the 
blue color of V2 would imply a reasonably promi- 
nent disk (the source of flickering). However, no 
evidence for flickering is seen as the residuals af- 
ter subtracting the sinusoidal model are consistent 
with the noise properties of the data. 

Besides being blue, variable objects, LMXBs 
also typically have weak emission lines in the op- 
tical and infrared. Shahbaz et al. (1996) show 
that the average H/3 equivalent width (EW) of 15 
LMXBs, is 1.9 A. Scaling by the average EW ratio 
between H/3 and Pa for the two CVs PQ Gem and 
EX Hya (using values from Dhillon et al. 1997; no 
published Pa EWs for LMXBs were found) we ex- 



pect an average Pa EW of 3.8 A. This is equivalent 
to F187N-F187W = -0.07, well below the la error 
for the faint end of our F187N vs F187N-F187W 
CMD (containing V2) of 0.11 mag. Therefore the 
Pa data is inconclusive in determining whether V2 
is an LMXB. 

An obvious constraint on the identify of V2 is 
its position relative to the X-ray source identified 
by Johnston et al. (1995). Our measured posi- 
tional offset of 5.5" is just outside the 5" error cir- 
cle quoted by Johnston et al. (1995). We hope to 
significantly reduce the size of this error circle with 
the use of our Cycle 1 Chandra observations. With 
these deep, 50ks observations the LMXB should 
have an error circle < 1" and should be distin- 
guishable from the large expected population of 
CVs and MSPs in this cluster, unless it is in deep 
quiescence. 

We may also compare the period and luminos- 
ity of V2 with UV or optically identified glob- 
ular cluster LMXBs. Three of the six globu- 
lar cluster LMXBs with optical counterparts have 
My >4.5, according to Deutsch, Margon and An- 
derson (2000). Figure 4 shows that this is be- 
yond our detection limits, for stars on the main 
sequence or the giant branch (for hotter objects 
our detection limits are even shallower). To make 
matters worse the Terzan 5 LMXB, being a tran- 
sient should be fainter than most other cluster 
LMXBs. We can predict My using the depen- 
dence of optical luminosity upon X-ray luminosity 
and period (van Paradijs and McClintock 1994) a 
relation that holds for both cluster and field sys- 
tems as shown by Deutsch et al. (2000). For a 
period of ~14 hrs, we estimate My = 4.2, beyond 
our detection limits. Finally, we note that 4 of the 
5 systems with optical counterparts and measured 
periods have period < 5.7 hours and three have 
periods < 0.73 hours. 

To summarize, although we find no individual 
argument gives compelling evidence against V2 
being the infrared counterpart of the Terzan 5 
LMXB, the accumulated weight of evidence argues 
against it. An alternative hypothesis is that V2 is 
an eclipsing BS. This is consistent with: (1) The 
position of V2 in the F110W vs F110W-F187W 
CMD (see figure 5). V2 has about the same lu- 
minosity as the brightest BS in 47 Tuc, and its 
luminosity may also be significantly enhanced by 
its binary companion. (2) The light curve ampli- 
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tudes of V2. Sinusoidal light curves with our mea- 
sured amplitudes are easy to reproduce using the 
light curve models of Wilson for W UMa systems 
(van Hamme, private communication). A growing 
number of binary BSs are being discovered in glob- 
ular clusters (eg Edmonds et al. 1996 and Kaluzny 
et al. 1998), although the fraction of BSs that are 
binaries is typically quite low. (3) The inferred 
orbital period of 14 hours is consistent with the 
periods of other binary BSs in globular clusters 
such as 47 Tuc (Edmonds et al. 1996). As with 
the LMXB hypothesis, the Pa data is inconclu- 
sive, since for a model BS with T e ff = 8000 K 
(at the top of our artificial BS sequence), F187N- 
F187W = 0.04, well within our 0.11 mag scatter 
for this line discriminant. 

Using our simulations we can ask how many 
eclipsing binaries we should expect to see with our 
data set. Since most eclipsing binaries discovered 
in clusters are found near the MSTO (eg Edmonds 
et al. 1996) and we have very little sensitivity to 
such variables, the simple answer is none. There- 
fore, a naive calculation may predict a very large 
number of eclipsing systems (and BSs) should be 
found in Tcrzan 5, if V2 is indeed an eclipsing BS. 
While clearly small number statistics apply since 
we have detected only one system, it is striking 
that evidence exists for a bright eclipsing BS at 
the edge of our detection limits. A relatively large 
number of BSs and binaries would have interest- 
ing consequences for the dynamics of this extreme 
globular cluster. For example, a high BS number 
may be caused by a high stellar interaction rate, 
which itself may be enhanced by a large binary 
fraction. A high interaction rate is expected for 
Terzan 5 because of its moderately sized core and 
extremely high stellar density (Cohn et al. 2000). 
Terzan 5 also may be the globular cluster with the 
largest population of MSPs in the galaxy (Fruchtcr 
& Goss 2000), further evidence for possible large 
interaction rates. 

It is useful to discuss the potential for detailed 
studies of other interaction products or binaries in 
Terzan 5. Our upcoming GO Chandra observa- 
tions of Terzan 5 should allow a deep search for 
CVs in this cluster, unless the LMXB is in out- 
burst. Unfortunately optical or IR identifications 
of such sources will be extremely difficult or impos- 
sible. We have already seen the problems with IR 
observations, where crowding prevented us from 



reaching even close to the MSTO, but with their 
smaller PSFs, Gemini and NGST will offer sub- 
stantial improvements over HST for studying this 
cluster. In optical passbands such as V, the PSF 
will be substantially reduced in size but the enor- 
mous reddening for this cluster implies that the 
MSTO will be found near V ~25.5, and therefore 
CVs at V > 28. Studies of the BS sequence will 
be difficult but may be feasible depending on the 
results of detailed crowding simulations. 
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-Noise limit (F187N) 




Fig. 1. — Plots of rms versus magnitude for the 
three different NIC2 filters used. The solid line 
shows the Poisson limits for aperture photometry. 
The symbols represent the variable VI ('1') and 
its near neighbors ('a'-'e'). 
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F110W 

Fig. 2. — Contamination level plotted against 
F110W. The same symbols as Figure 1 are used 
except we also include the variable V2 ('2'). 




Fig. 3. — The F110W image for the first orbit, 
including the positions of VI and V2 and the po- 
sition and error circle of the X-ray source from 
Johnston et al. (1995). The close-up figures in- 
clude solid circles for the variables and dotted cir- 
cles for the near neighbors of VI ('b'-'e'; 'a' is not 
circled and is found 2 pixels above VI) and V2. 
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Fig. 4. — The percentage of detected simulations 
(with 0.3 mag amplitude) as a function of radial 
distance from the center of the cluster. A range 
of absolute magnitudes were used for the simula- 
tions. The drop in detection rate beyond a radial 
distance of ~ 10" is because of the decreased frac- 
tion of annulus that lies in the NIC2 field. 
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Fig. 5.— F110W vs F110W-F187W CMD from 
Cohn et al. (2000). The numbered and labeled 
stars discussed in the text are shown, along with a 
selection of blue stars ('*')• Isochrones from Berg- 
busch and Vandenberg (1992) are shown by the 
solid line, plus HB and BS sequences. The jagged 
appearance of the isochrones is because they have 
been fed through Kurucz atmospheres, quantified 
in T e ff and log g, to calculate the NICMOS mag- 
nitudes. 
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Fig. 6.— Time series for VI plus 5 near neigh- 
bors. The unfilled squares show the F110W data 
(averaged over each orbit), the filled squares the 
F187W data (orbital average) and the solid lines 
the F187N data. For VI the zero-points have been 
adjusted so that the F187W and F187N filters 
have the same magnitude for orbit 3 as F110W. 
For the neighboring stars the zero-points have 
been adjusted so that the average magnitude of 
the three filters over all of the orbits are equal. A 
further correction was applied to offset the time 
series from each other. 



Fig. 7. — Time series for the RR Lyrae variable 
(VI). The unfilled squares show the F110W data, 
the filled squares the F187W data and the solid 
lines the F187N data. The fitted J, H and K 
templates are shown by the dotted, dashed and 
dot-dashed lines respectively. 



15 



(a) 



o 



12 - 



14 



16 - 



20 



(b) 



10 



12 



co 

E 14 



16 - 



. I. 



2 



ffi- 



+ + 



1 2 3 

F110W-F187W 



-0.4 -0.2 



0.2 



0.4 



F187N-F187W 
Fig. 8. — CMDs from combining the time se- 
ries photometry. Figure 8(a) shows the F110W 
vs F110W-F187W CMD and Figure 8(b) shows 
the F187W vs F187N-F187W CMD. Only stars 
with CL < 1.5 in F187W have been plotted here. 
The numbered and labeled stars from Figure 5 are 
shown along with extra blue stars from our pho- 
tometry ('+'). 
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Fig. 9. — A similar plot to Figure 1 except that 
only stars with CL (F187W) < 1.5 have been plot- 
ted. The variable V2 has been added to this figure, 
along with a symbol 'R' showing the residual after 
subtracting a sinusoidal fit to the V2 light curve. 
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Fig. 10. — Light curves for V2. The upper panel 
shows the aperture photometry results (including 
a sinusoidal model fit to the data), the next panel 
the noiser PSF-fitting results, then the residuals 
after subtracting the model. The other panels 
shows the light curves of 5 nearby stars, with the 
required F110W offset shown in parentheses. 




Fig. 11.— A finding chart for the NIC1 MSP field. 
The nominal MSP position is shown (M'), along 
with HBs ('H') and red giant stars ('R') discussed 
in the text. A 1" circle is plotted around the nom- 
inal MSP position. 
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Fig. 12. — The plot of rms versus magnitude for 
the NIC1 field containing the MSP. Stars within 
1" of the nominal MSP position are shown, along 
with HB stars and red giant stars discussed in 
the text. The estimated position of the turnoff 
is shown. 
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